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ABSTRACT 

Context. The progenitors of many Type II supernovae have been observationally identified but the search for Type Ibc supernova (SN 
Ibc) progenitors has thus far been unsuccessful, despite the expectation that they are luminous Wolf-Rayet (WR) stars. 
Aims. We investigate how the evolution of massive helium stars affects their visual appearances, and discuss the implications for the 
detectability of SN Ibc progenitors. 

Methods. Evolutionary models of massive helium stars are analysed and their properties compared to Galactic WR stars. 

Results. Massive WR stars that rapidly lose their helium envelopes through stellar-wind mass-loss end their lives when their effective 

temperatures - related to their hydrostatic surfaces - exceed about 150kK. At their pre-supernova stage, their surface properties 

resemble those of hot Galactic WR stars of WO sub-type. These are visually faint with narrow-band visual magnitudes M r = -1.5 

2.5, despite their high bolometric luminosities (log L/L Q = 5.6 ■ • • 5.7), compared to the bulk of Galactic WR stars (M,. < -4). In 
contrast, relatively low-mass helium stars that retain a thick helium envelope appear fairly bright in optical bands, depending on the 
final masses and the history of the envelope expansion during the late evolutionary stages. 

Conclusions. We conclude that SNe Ibc observations have so far not provided strong constraints on progenitor bolometric luminosities 
and masses, even with the deepest searches. We also argue that Ic progenitors are more challenging to identify than lb progenitors in 
any optical images. 

Key words. Stars: evolution - Stars: Wolf-Rayet - Stars: binary - supernovae: general 



1. Introduction 

Direct detections of the progenitor stars of supernovae discov- 
ered in the nearby Universe provide one of the most stringent 
constraints on stellar evolution theory. Since the unambiguous 
identification of the progenitor of the supernova 1987A, ob- 
servers have identified progenitor stars of numerous Type IIP 
supernovae and severa l Type lib an d Iln supernovae in pre- 
supernova images (see ISmarttl 120091 for a review). However, 
despite their relevance to Galactic chemical evolution, the pro- 
genitor stars of supernovae (SNe) Ibc remain as yet elusive. 

For decades, there had been a widely held belief that SNe 
Ibc progenitors are massive Wolf-Rayet (WR) stars, formed ei- 
ther through stellar-wind mass-loss or Roche-lobe overflow in 
a binary system. Also the identification of broad SN Ic fea- 
tures in the vicinity of several long GR Bs dGalama et al. Ill998t 
iHjorth et al.ll2003l: IStaneket al.l 120031) for which th e progen- 
itors were also suggested to be massive WR stars dWooslev I 
fl993h added further evidence to this assertion. An alterna- 
tive to t he massive WR scenario i s that of lower mass bina- 
ries (e.g.lPods adlowski et al. lll992r.lWellstein & Langerlfl999t 



supernova stage dMevnet & Maederll2003l:ICrockett et al. 120071: 



lEldridge et al. ll2008tlYoon et al. 120101) 



Although many previous searches for SNe Ibc progenitors 
were hampered by insufficient detection limits, WR stars with 
log L/L Q > 5.3 could have been i dentified as the progeni 



tors of the Type lb supernova 2000ds I Maund & Smartt Il2005h . 
and the Type Ic supernovae 2004gt dMaund et al. 1 120051) and 
2002ap dCrockett et al. 1120071) . On the other hand, the most re- 
cent stellar evolution models predict that WR stars originating 
from massive single stars with initial masses higher than about 
30 M Q have bolometric luminosities of log L/L Q > 5.4 at the pre- 



iGeorgv et al. (12012 ). New detection limits, e.g., on SNe 2002ap, 
have mostly been interpreted as evidence against the massive 
WR scenario, but in favour of progenitors of lower mass binaries 
inste ad and largely dismissing the option of massive WR stars 
(e.g. JCrockett et al. 112007b ISmartt ll2009HYoon et al. Il2010h . 

However, in these interpretations the crucial final stages of 
massive-star evolution are not properly accounted for, as we ar- 
gue in the following. In this paper, we present the seemingly 
counter-intuitive idea that the more massive - and bolometrically 
more luminous - single WR progenitors are more challenging to 
detect than low-mass binaries. 



2. Massive Wolf-Rayet stars: Models and 
observations 

One of the most crucial questions regarding the observational 
identification of massive WR stars as SN Ibc progenitors is this: 
which bolometric correction (BC) should be applied to WR stars 
at the pre-supernova stage? Because the optical luminosities of 
WR stars are critically affected by optically thick WR winds, the 
simple assumption of black-body ra diation should no t be applied 
eve n for a crude approximati on (cf .lCrow ther 20 071). 

lMaund& Smartt I (120051) and iMaund et all d2005l) calcu- 
lated BCs f or WR stars using the Potsdam grid of synthetic 
WR spectra (lHamann & Grafener 1 l2004t iGrafener et al. ||2002|) 
and obtained the detection limits for different SNe Ibc progeni- 
tors in terms of bolometric luminosities. These detection limits 
were the n compared to the ev olutionary tracks of massive stars 
given by iMevnet et al. I (1 19941) to constrain the SN Ibc progeni- 
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tor masses. The caveat here is that M evnet et al. I (1 19941) applied 
a correction to the radii of their WR star models to take into ac- 
count the radial extension of th e sta r due to an opt ically thick 
wind, following lLangerl (1 19891) and lMaederl dl990l) . This cor- 
rection is particularly large as they employed high WR mass-loss 
rate s for which the effect of win d clumping was not considered 
(cf. lHamann & Koesterkelll998l) . As a consequence, the evolu- 
tionary tracks of their WR star models have surface temperatures 
of log T - 4.8, and Maund et al. applied their BCs accordingly 
only up to the quoted value. However, the stellar temperatures 
T+ of the Potsdam models do not refer to the values at the pho- 
tosphere within the extended optically thick wind, but to those 
at large optical depths close to the hydrostatic stellar surface. 
This inconsistency introduced a large underestimate of the BCs 
applied to WR stars at the pre-supernova stage, because stel- 
lar evolution models predict that massive WR stars with masses 
M > 10 M Q end their lives with hydrostatic surface tempera- 
tures well above log T+ — 5.0 (s ee discussion below). On the 
other hand, ICrockett et al. I (120071) simply employed a constant 
BC of -4.5 for the ent i re tem perature range of WR stars, follow - 
ing lSmith & Maederl d 19891) . More recent analyses of WR stars, 
however, seriously question the validity of this approach, as dis- 

cus sed below. 

lHamann et aTl d2006l) and ISander et al71 d2012l) (hereafter, 
the Potsdam group) provide a homogeneous set of data on 
the properties of Galactic WR stars of different spectral types 
(WNL, WNE, WC, and WO). Fig.Q]shows the absolute narrow- 
band visual magnitudes M,Q of the Galactic WR stars with 
known distances, as given by the Potsdam group. We note the 
strong temperature dependence of M v : WR stars tend to be visu- 
ally dimmer at higher temperatures. The corresponding BCs are 
presented in Fig. [2] for which the best linear fit is given byQ 



BC = 22.053 - 5.306 log . 



(1) 



To compute BC = M bo \ - M v , we adopt M v - M v = 0.75 to cor- 
rect for the effect of stro ng emission lines on the visual broad- 
band ma gnitude M v (cf. ICrockett et al. 112007b . The BCs calcu- 
lated by iMaund et al. I (12005 ) using the Potsdam WR grid are 
plotted in the same figure for comparison. Here again, the figure 
reveals that BCs for WR stars are a sensitive function of tem- 
perature, and that a more negative BC should be applied at a 
higher temperature. The synthetic spectra from the Potsdam WR 
models agree reasonably well with observations, confirming the 
temperature dependence of BCs. 

To determine the implications of this finding for SN Ibc pro- 
genitors, we show a Hertzsprung-Russell (HR) diagram of evo- 
lutionary tracks of single massive helium stars of solar metallic - 
ity in Fig. [3] Our models were calc ulated with the hydr odynamic 
stellar-evolution code des cribed in Yoon et al I (20101) using the 
WR mass-loss rates from iNugis & Lamers I (120001) . The figure 
shows that helium stars with initial masses higher than about 



1 Throughout this paper, M v represents an absolute magnitude in the 
narrow visual-band (3900 <A< 4500). 

2 In stellar evolution models, the hydrostatic surface temperature 
corresponds to the effective temperature at the surface boundary for 
which radial extension of the star due to optically thick WR winds is not 
considered. The empirical stellar temperatures T+ in the Potsdam group 
data we re estimated at large optical depth (t Ross = 20) in the WR wind 
models dSander et al7ll2012h and thus nearly resemble the hydrostatic 
values. However, as the hydrostatic layers are not directly observable, 
the empirical rely on the adopted wind structure. For the stars with 
the strongest winds, we estimate a possible systematic uncertainty of up 
to 0. 1 in log caused by this effect. 
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Fig. 1. Absolute visual magnitudes of Galactic Wolf-Rayet stars 
with known distances of different spectral types (WNL, WNE, 
WC and WO) as a function of surfac e te mperatures. The dat a 
were taken from lHamann et al~ ] (120061) and lSander et"al~l d2012l) . 
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Fig. 2. Bolometric corrections (BCs) for Galactic WR stars, cal- 
culated from the data given in Fig. [1] as a function of surface 
temperatures. The thick dashed line gives t he best fit to the data . 
The thin lines denote the BCs calculated bv lMaund et al. ( 20051) 
from t he Potsdam WR-star synthetic spectra of Grafener et al. I 
d2002l) . for the minimum (dotted line), median (dashed line), and 
maximum (dot-dashed line) transformed radii R, (see Maund et 
al. for details). 



15 M become gradually hotter. This is because they progres- 
sively lose their helium-rich envelopes in a stellar-wind, expos- 
ing the carbon-oxygen cores that gradually contract. The cor- 
responding final masses are higher than about 9.0 M Q and the 
surface temperatures increase to log 7V 5.3 at the end of 
their evolution. We note that the recent single star-models of 
the Geneva group predict similar surface temper atures for WR 
stars of log L/L Q > 5 .3 at the pre-supernova stage dGeorgy et al. I 
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120121) . The question to be addressed now is which bolometric 
correction should be applied to such WR stars with log L* > 5.3. 

Fig. [3] shows the well-known discrepancy between the pre- 
dicted surface temperatures of WR stars and the observation- 
ally derived values: the helium-star models gener ally give higher 
tempe ratures than those of observed WR stars. IGrafener et al. I 
d2012l) suggested that the observationally implied large radii of 
many Galactic WR stars are explained by envelope inflation 
caused by the iron opacity peak at logL =* 5.18 and the ef- 
fect of density inhomogeneities (clumping) in the inflated lay- 
ers. These inflated envelopes usually contain a str ong density 
inversion, as is often obse rved in stellar models (llshii et al. I 
1 19991: IPetrovic et al71l2006l) . This density inversion is also ob- 
served in our helium-star models presented in Fig. [3] but the 
degree of infla ti on in our models is much weaker than in the 
IGrafener et al. I d2012l) models because we do not consider the 
effect of clumping. Interestingly, given the key role of the iron 
peak, this inflation may no longer occur when the surface tem- 
perature of a star approaches or becomes hotter than the iron 
peak temperature (log T - 5.18). Our helium-star models in- 
deed show that the density inversion in the outermost layers, 
a sign of inflation caused by the iron peak, disappears when 
log T+ > 5.07. This implies that, although massive WR stars 
have inflated envelopes for most of their lifetimes, the infla- 
tion disappears as they lose most of their helium envelopes and 
their c ores become more compact during the late evolutionary 
stages. Grafe ner & Hamannl d2005l) showed that these compact 
stars form hot WR-type winds driven by the radiative force on 
the iron peak opacities. This may explain why the envelopes of 
these stars do not inflate and the positions of early-type WC and 
WO stars agree well with the model predictions, while WR stars 
of later spectral subtypes have much cooler temperatures. 

Although we cannot directly apply the BCs of Fig. [2] to the 
stellar evolution models, we can draw the following conclusion: 
massive WR stars experiencing strong mass-loss end their lives 
at their hottest point (log L+ = 5.3). Therefore, SN Ibc progen- 
itors with final masses higher than about 9 M e resemble the 

WO stars with log T+ =* 5.2 and M v = -1.5 2.5 (Fig. \J}, 

rather than the majority of observed WR stars (log L* < 5 and 
M v < -4). They are relatively faint in optical bands compared to 
the majority of observed Galactic WR stars, and we cannot ex- 
clude massive WR stars of log L/L > 5.4 as the progenitors of 
SNe 2000d s, 2004gt and 2002ap, in contrast to the p r evious con- 
clusions oflMaund & Smarttl (l2005h . iMaund et all (120051) . and 
ICrockett et al. I (120071) . 

On the other hand, the helium envelopes in less-massive 
helium stars (M = 8 and 10 M Q Fig. [3} are not greatly 
stripped off and carbon is not enriched at the surface. These 
stars therefore rapidly expand during the late evolutionary 
stages b eyond core helium-burn i ng be cause of the mirror ef- 
fect (cf. iKippenhahn & Weigert I [19901) . The surface tempera- 
ture decreases from log L* 5.2 at core helium-exhaustion to 
log T+ =* 4.9 at core oxygen-burning, accordingly. Although 
the bolometric luminosities of these SN Ibc progenitor mod- 
els (log L/L « 5.2) are lower than those of more massive 
ones (log L/L Q > 5.4), much lower temperatures at the pre- 
supernova stage imply that they should appear more luminous 
in optical bands. Given that WR stars with log L* < 5.0 and 
logL/Lo < 5.3 in the Potsdam sample have M v - -4 - ■ • - 5, 
similar M v s are expected for these SN Ibc progenitors. 
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Fig. 3. Evolutionary tracks of pure helium-star models at solar 
metallicity in the HR diagram, for six different initial masses (8, 
10, 15, 20, 25, and 30 M Q ), as indicated by the labels. The num- 
bers in the parentheses denote the corresponding final masses. 
The positions where the surface mass-fraction of carbon is larger 
than 0.2 are marked by thick grey lines. The star symbol marks 
the end point of the evolution, which is the beginning of core 
oxygen-burning that occurs a few years before core collapse, for 
each evolutionary sequence, while the tilted rectangle denotes 
the position at 1000 yrs before the supernova explosion. The 
positions of G a lactic WR stars with known dis tances given by 
lHamann etaTI (120061) and ISander et all (1201 2l) are marked by 
open squares (WNL), filled squares (WNE), filled circles (WC), 
and triangles (WO). 



3. Relatively low-mass helium stars in binary 
systems 

Galactic WR stars have masses higher than about 8 M e and 
many of them are thought to have originated from massive single 
stars through stellar-wind mass-loss. In binary systems, lower- 
mass heli um stars can be pr oduced via mass transfer. As dis- 
cussed bv lYoon et al~l d2010h . such relatively low-mass helium 
stars can retain a large amount of helium in the envelope and ex- 
pand to a large radius during the late evolutionary stages beyond 
core helium exhaustion. The final radii of SN Ibc progenitors are 
hence systematically larger for lower masses. 

Fig. |4] shows the evolutionary tracks of several helium star s 
produced in close binary systems taken from lYoon et al. I (120 1 Oh . 
Some of the helium stars in this figure (YWL9, YWL14 and 
YWL21) end their lives while they are filling their Roche- 
lobes during Case ABB or BB mass transfer (i.e., mass trans- 
fer from the He s tar during/after th e carbon-oxygen core con- 
traction phase; see I Yoon et al. II20T0I for more details). For com- 
parison, the evolutionary track of a single, pure helium star of 
M = 3 M Q is also presented in the figure. Since these relatively 
low-mass helium stars may not have optically thick winds, the 
black-body assumption is adopted as a rough approximation to 
calculate the corresponding visual magnitudes, which are shown 
in Fig. [4] This figure shows that helium stars of M — 3M -5M Q 
are visually very faint (M v > 1) on the helium main-sequence. 
However, rapid expansion of the helium envelope during and af- 
ter the carbon-burning phase suddenly makes these stars lumi- 
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Fig. 4. Upper panel: Evolutionary tracks of helium stars pro- 
duced in binary systems through Case AB or B mass transfer 
phase, and of a single helium star of 3 M G (black solid line), at 
solar metallicity. The star symbol marks the end point of each 
evolutionary track, which is the core carbon burning phase for 
3 M Q star and the oxygen burning phase for the others. The 
helium -star models in binary systems are taken from lYoon et al. I 
d2010l) . and the evolutionary tracks are plotted only from the 
onset of core helium burning. The sequence number given for 
each track (except for the 3 M Q single helium star) denot es the 
model sequence number in Table 1 of I Yoon et al. (l2Q10h . The 



final mass (Mj) for each model sequence is given in the paren- 
thesis. Lower panel: Corresponding absolute narrow-band visual 
magnitudes calculated with the assumption of black-body radia- 
tion. 



nous in the visual band shortly before their deaths. In particu- 
lar, the surface temperature of 3 M Q helium star can decrease to 
log T+ =3.8 and it may be particularly luminous in the R band 
(M R - -5.5). 

The important implication of this result is that some rela- 
tively low-mass SN Ibc progenitors (~ 3 M Q ) are more lumi- 
nous in either the visual band (M v -4) or other optical bands 
than much more massive SN Ibc progenitors with final masses 



higher than about 9 M Q , which should be as dim as M v ^ -1.5 
as discussed above (see also Fig.Q]). 

We note that not all SNe Ibc progenitors with final masses of 
M < 5 M Q may be optically bright. Because they are likely to 
undergo mass transfer in a close binary system during the final 
evolutionary stages owing to the rapidly expanding envelopes, it 
is possible that some of them may lose most of their helium en- 
velopes by the time of su pernova explosion (if the orbital separa- 



tion is sufficiently short: iPols & Dewi 



Ivanova et al. 1 120031: lYoon et aTT 201C ). These relatively low 



2002: lDewi etaLl [2002: 



mass, helium-deficie nt stars may constitute the low-mass class 
of SN Ic progenitors dPols & Dewi ll2002l:lYoon et al. 1201 d) and 
may explain the fast light-curve of SN 19941 ( Nomoto et al. I 
|1994|) . Such a SN Ic progenitor of the low-mass class has a sur- 
face temperature of log > 5.0 and a very faint optical lu- 
minosity (M v > 1.0). On the other hand, helium stars of about 
2 M Q - 3 M Q can also be produced by a merging of a white dwarf 
and a helium star, in which case their surface temperatures may 
decrease below 10 4 K at the pre-supernova stage, which would 
appear to be optically bright single stars. Therefore, fairly di- 
verse properties are expected for SN Ibc progenitor systems of 
this mass range. 

4. Discussion 

Our study refutes the widespread belief that more massive pro- 
genitors of SNe Ibc are more easily identifiable than less massive 
ones. In contrast, our analysis implies that massive SN Ibc pro- 
genitors with final masses M > 9 M and bolometric luminosi- 
ties log L/Lq > 5.4 may be optically faint compared to relatively 
low-mass SN Ibc progenitors, depending on their mass-loss his- 
tory. 

The key factor that determines the brightness of a SN Ibc 
progenitor in optical bands is the mass of the helium envelope. If 
most of the helium envelope is lost via either stellar winds or bi- 
nary interactions, the SN Ibc progenitor will be very hot, result- 
ing in a low optical luminosity at the pre-supernova stage. This 
may be the case for the most massive WR stars as discussed in 
Sect. [2] If a fairly thick helium envelope can be retained in a SN 
Ibc progenitor, the expansion of the envelope that occurs shortly 
before the SN explosion makes the progenitor star fairly bright 
in optical bands. This case is mostly relevant for relatively low- 
mass helium stars produced in close binary systems (Sect. [3). 
Therefore, we suggest that SN lb progenitors would be optically 
brighter - hence easier to detect - than SN Ic progenitors, in 
general. 

We summarize in Table Q] our rough predictions for the vi- 
sual magnitudes of SNe Ibc progenitors for different possible 
cases. A secure identification of a SN Ic progenitor would be 
challenging. A detection limit of M v - -1.5 for SNe Ic of the 
high-mass class (i.e., Mj > 9 M ) and a much better detection 
limit is needed fo r SNe Ic of the l ow-m ass class. Therefore, it is 
unsurprising that lCrockett et al.1 d2007l) could not find the pro- 
genitor star of SN Ic 2002ap even with the best detection limit 
achieved so far (i.e., M B > -4.2 and M R > -5.1). Although SNe 
lb progenitors are predicted to be systematically more luminous 
than SNe Ic progenitors in optical bands, their visual luminosi- 
ties would still be limited to M v - -5 and a wide range of optical 
magnitudes are expected. 

Given that many SNe Ibc are expected to occur in binary sys- 
tems, another factor that should be considered in the search for 
SNe Ibc progenitors is the companion star. Even if a SN Ibc pro- 
genitor is optically faint, the companion star may be detectable. 
The binary star models indeed predict that a significant fraction 
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Table 1. Predicted visual magnitudes at the pre-supernova stage for different final masses of helium stars. 



Final Mass 


SN Type a 


Helium envelope 


Single/Binary 


My 




M f > 9 M G 


Type Ic 


Stripped-off 


Single or Binary 


My = -1.5- 


••-2.5 


6 < M f < 8 M B 


Type lb 


Retained 


Single or Binary 


My = -4 • • ■ 


-5" 


2 < Mf < 5 M Q c 


Type lb 


Retained 


Binary or merger remnant 


My =-!■■■ 


- 4 (M R = -2.0 5.5) 


2 < Mf < 5 M c 


Type Ic 


Stripped-off 


Binary 


My > 1 





Notes. <n> SN Ic is assumed if the helium envelope is stripped off from the progenitor (M ^ $ 0.5 M Q ) and SN lb otherwise. However, this may 
also depend on the history of chemical mixing during the supernova explosion as shown by D essart et al. I d2012f) . This estimate is based on the 
assumption that the helium stars of this mass range still have optically thick WR winds. If not, M v is larger (M r > -2). w The lower limit of 2 M e 
is an arbitrary choice, but SNe Ibc with M < 2 M Q are likely to be very faint and difficult to discover. 



of SNe Ibc progenitors have an O-type companion. However, 
not all of them are bright enough to be identified with the de- 
tection limits achieved in the previous observations. Our prelim- 
inary calculation indicates that only a small fraction (less than 
about 10%) of SNe Ibc progenitors have a luminous O star com- 
panion with M v < -4. 

5. Conclusion 

Our study has shown that the masses of SNe Ibc progenitor stars 
are not linearly correlated with their optical brightness and that 
the evolution of the surface properties of massive helium stars 
during their final evolutionary stages should be carefully inves- 
tigated. The non-detection of a SN Ibc progenitor even with a 
good detection limit does not necessarily imply that its progen- 
itor is a relatively low-mass helium star, and vice versa. This 
should be properly taken into account in future observational ef- 
forts to directly identify SNe Ibc progenitors. 
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